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ABSTRACT 

Near surface measurements of turbulent temperature and 
velocity fluctuations in the ocean are obscured by relatively 
much larger ordered fluctuations due to surface gravity waves. 
Analysis of measured velocity signals in the near surface 
region indicated that the ratio of turbulent to wave-induced 
kinetic energy was consistently less than 10" . To a first 
approximation, wave-induced fluctuations should not contribute 
to turbulent transport. However, small instrument errors due 
to misalignment, phase response, or directional response may 
cause significant wave-induced error in direct flux measure- 
ments. The effect on turbulent flux measurements of misalign- 
ment of the velocity sensors was examined. It was found that 
misalignment resulted in apparent phase shifts in the measured 
wave-induced signals, causing erroneous contributions to the 
calculated flux. For alignment errors of 3.6°, the error in 
calculated momentum and heat fluxes could have been up to 
500% abd 110%, respectively. Momentum flux measurements made 
in deep water (with respect to the length of surface waves) 
were found to be less sensitive to alignment errors. The 
errors in this case were functions of the alignment of the 
flowmeter axes with the predominant wave direction. 
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I, INTRODUCTION 



Turbulent processes in the upper ocean are not yet fully 
understood, Phillips (1966), Kitaigorodskii (1971) and others 
have pointed out the need for near surface measurements of 
turbulence and turbulent transfer. Large coordinated 
experiments such as BOMEX have been designed to measure 
turbulence in the Marine Boundary Layer of the atmosphere 
(Holland, 1972), but complementary measurements in the near 
surface layer of the ocean have been scarce. The recently 
completed GATE should greatly add to such data. Several 
projects (NORPAX, JASIN, GARP) are currently being planned 
to measure turbulence in what is generally termed the mixed 
layer of the ocean. 

The dynamics of the near surface layer (in the absence 
of tropical storms) are dominated by surface gravity waves. 

The ordered motion of the waves is usually not of interest 
when making turbulence measurements because they do not 
contribute, at least to a first approximation, to the vertical 
transfer of heat and momentum through the water column, As 
a result, near surface turbulence measurements seek to 
resolve what are essentially second order quasi-random 
motions against a background of first order coherent noise. 

The Naval Postgraduate School has conducted a series of 
turbulence measurements using a fixed array of sensors 
mounted on the Naval Undersea Center oceanographic research 
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tower off San Diego (Hagen, 1974), An attempt was made to 
measure turbulent heat and momentum fluxes from simultaneously 
recorded time series of orthogonal velocity components and 
temperature. It was found that wave-induced fluctuations 
may have caused significant errors in the measurements. In 
view of the current emphasis on near surface turbulence 
measurements, an analysis of wave-induced measurement errors 
was considered appropriate. Only through such analyses can 
realistic experiments be designed. 

The analysis of wave-induced errors in turbulence measure- 
ment is carried out first, followed by a description of the 
experiment and the results obtained. 

II . WAVES, TURBULENCE AND NOISE 

In the upper layer of the ocean, turbulence is superim- 
posed on ordered motions related to surface gravity waves. 

The contribution of the waves to a measured signal may be 
quite large. Data taken from the NPS experiments were 
examined as to the relative magnitudes of wave-induced 
signals, turbulence and noise. The data were presumed repre- 
sentative of open ocean as well as coastal near surface 
environments , 

Fluctuations of horizontal and vertical velocity (u, w) 
and temperature, ©, can be thought of as composed of wave- 
induced and turbulent parts: 
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u = U + u f 

w = W + w* 

0 = 0 + e* (i) 

where upper case letters stand for wave-induced fluctuations 
and primed letters for turbulent fluctuations. With certain 
assumptions it is possible to separate the wave-induced and 
turbulent parts of equation (.3.) and estimate their relative 
magnitudes. For moderate wave conditions, it is reasonable 
to assume that the wave-induced spectral velocity components 
can be calculated from the wave spectrum using linear wave 
theory. The success with which wave velocities are 
described by linear wave theory can be gauged by applying 
the appropriate linear transfer function to the surface 
displacement spectrum and studying how well the calculated 
and measured velocity spectra compare. Figure (1) shows the 
measured vertical velocity spectrum at 6.6 meters depth and 
that calculated from the surface displacement spectrum using 
linear wave theory. Figure 2 shows the coherence and phase 
functions obtained from the cross spectrum between measured 
wave and velocity records. The maximum coherency squared 
(hereafter referred to simply as the coherence) occurs at a 
period of 5.9 seconds and is equal to 0.9948. The phase 
function shows the waves leading the vertical velocities by 
approximately 90° throughout the surface wave band. Figures 
1 and 2 clearly demonstrate the applicability of linear wave 
theory to wave-induced motions. 
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Figure 1, Comparison between measured vertical velocity 
spectrum and the vertical velocity spectrum 
calculated from the waves using linear theory. 
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Figure 2, 



Coherence and phase functions calculated from the 
cross spectrum between surface displacement and 
vertical velocity. 



It is assumed that wave-induced and turbulent motions 
are independent (Phillips, 1966) (_Seitz, 1971), The inter- 
actions between wave motions and turbulence contribute to 
the motion only at the third order (Kinsman, 1965), There- 
fore, the independence assumption is justified except for 
very steep waves. 

Using the two preceding assumptions, the measured velocity 
spectrum can be separated into wave-induced and turbulent 
parts using the experimentally determined coherence between 
the wave displacement and velocity signals (Benilov and 
Filyushkin, 1970, see Appendix A), The turbulence spectrum 
is then given by 

S ,(f) = S (f)(l - Y 2 ) (2) 

vr 7 w y 1 nw y v y 

and the wave-induced spectrum by 

s ff (f) = S w (f) - S w ,(f>. (3) 

The variance of velocity due to wave-induced or turbulent 
motions is obtained by summing their respective spectral 
components. Table 1 lists experimentally determined velocity 
variances from the NPS data. The horizontal wave velocities 
are larger than the vertical velocities due to the shallow 
water environment. The magnitude of horizontal "turbulent” 
variance may not all be attributable to turbulence due to 
the effect of directionality of horizontal wave motions on 
the coherence (Yefimov and Khristoforov , 1971), Table 1 
indicates that the intensity of turbulent motions decreases 
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with depth roughly as the intensity of the wave motions. The 
ratio of turbulent to wave-induced energy is seen to be 
consistently less than 10^. 

An inseparable part of the turbulent variance is that due 

to instrument noise. The advertised broadband noise level 

for the flowmeters used (Marsh-McBirney Model 711) was 1.0 
2 -2 

cm sec , The noise level was checked by recording the 
instrument output while the flowmeters were at rest in still 
water and was visually found to agree with that advertised. 
Additional noise can be introduced in the form of turbulence 
created by the obstruction of the flow field by the sensor. 

At the small velocities measured this effect can probably be 
ignored (NOIC, 1974). 

Making the additional assumption that temperature can be 
treated as a passive scalar, the assumptions made with regard 
to velocity fluctuations should apply equally well to the 
wave-induced and turbulent temperature fluctuations of 
equation (1). The wave-induced temperature fluctuations can 
then be written in terms of the vertical displacement of 
some local temperature gradient relative to the point of 
measurement (Thornton and Boston, 1974): 



99o = _W 36o 
3z ia 3z ' 



(4) 



where Z is the wave-induced vertical displacement and 
30o/ 3z is the temperature gradient. Equation (4) implies 
that the coherence between 0 and the surface displacement 
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is 1,0 and that the method of Appendix A can be used to 
separate wave-induced and turbulent temperature components. 
Table 2 shows the wave-induced and turbulent temperature 
variances calculated from the NPS data, 

III, APPLICABILITY OF TAYLORS HYPOTHESIS 

An associated problem of turbulence measurements in 
the near surface layer, although not directly tied to flux 
measurements, deals with the applicability of Taylor's 
hypothesis. Experimental turbulence data is usually in the 
form of a time history of turbulent fluctuations measured 
at a point or with a towed sensor. The spectra calculated 
from such measurements are then defined in frequency space. 
Unfortunately the theoretical treatment of turbulence is 
usually confined to wave number space, so that a trans- 
formation is employed. Turbulent frequencies are not uniquely 
related to their wavenumber, as is the case for wave motions 
for which a dispersion relationship exists. Instead, use 
is made of Taylor's hypothesis, or the frozen turbulence 
approximation. Taylor's hypothesis states that if a sensor 
is moved rapidly enough through a turbulent field, the time 
dependent changes in the field are negligible, with the 
result that it is the spatial variations that are being 
measured, A transformation from time to space can therefore 
be made by substituting x/U for t, where U is the velocity of 
the sensor (Tennekes and Lumley, 1972), For fixed point 
measurements the transformation can still be made, providing 
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that the turbulent field is being swept by the sensor fast 
enough, as by a mean wind or current. The criterion for the 
applicability of Taylor's hypothesis is that the relative 
velocity of the sensor with respect to the fluid be much 
larger than the magnitude of turbulent velocity fluctuations: 

u'/ U << 1, 

Meteorologist never seem to have trouble making the 
transformation because mean wind speeds are on the order of 
meters per second while small scale turbulent fluctuations 
are on the order of centimeters per second. Mean flow 
velocities in the ocean are usually measured in centimeters 
per second, as are turbulent fluctuations. In the ocean it 
is necessary to verify whether or not Taylor's hypothesis 
is valid, or to tow the sensor through the water at a speed 
considered sufficient to satisfy the approximation. The 
latter method has been used most often to date. 

Seitz (1971) and Bowden (1963) concluded Taylor's 
hypothesis was valid in a tidal estuary for frequencies down 
to 0,25 Hz when the mean convective velocity was on the order 
of 25 cm/sec"' 1 '. The NPS experiment did not include a good 
measurement of mean velocity, but turbulent fluctuations were 
on the order of 0.3 cm sec”'*' per spectral estimate (Af = 
0.0067 Hz). A mean velocity of only 3.0 cm sec”'*' could 
conceivably be sufficient to satisfy Taylor's hypothesis in 
this case. 
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High mean velocities are seldom encountered in the open 
ocean. Even where surface currents may be large enough, the 
measuring platform usually moves with the current so that 
relative velocities might be negligible. In such cases, it 
is doubtful whether Taylor^s hypothesis can be applied. 



IV . THEORY OF TURBULENT FLUX MEASUREMENT 



If the Navier-Stokes equation is rewritten using mean 
and fluctuating velocities and then time averaged, one of 
the resulting terms is the Reynold's stress tensor (Phillips, 
1966) : 



- P u.u . 
i 3 



i, 3 = 1, 2, 3 



(5) 



were u^ are velocity fluctuations and the overbar denotes a 
time average. The off diagonal terms tif j ) , represent a 
turbulent transport of momentum and are covariances of 
orthogonal velocity fluctuation components. If the motion 
can be considered horizontally homogeneous, the downstream 
component of vertical momentum flux is given, in conventional 
notation, by 

T 

u(t)w(t) dt, (6) 

An analogous expression for vertical heat flux can be 
derived from thermodynamic relations and is 

T 





0(t)w(t) dt , 



(7) 
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where 0 represents a fluctuating temperature. The covariance 
between two fluctuating quantities can also be expressed 
(Bendat and Piersol, 1971) 



where Co^Cf) is the real (.in phase) part of the u-w cross- 
spectrum, The spectral representation is the most popular 
because turbulence theory has envolved in the spectral domain. 
Experimental determination of heat and momentum flux can be 
made by recording fluctuating velocity components u and w, 
and temperature fluctuations, 0 , at a point and calculating 
the covariance or cospectrum as above. 

Alternate methods for estimating turbulent heat flux 
are the dissipation method used by Williams and Gibson 
(1974)^ and the profile method. The latter is essentially 
an integration over time of an assumed equation of conser- 
vation for heat (Voskanyan, et, al . , 1969) (Kaiser and 
Williams, 1974). Both these methods have the advantage that 
only temperature measurements are necessary, while the eddy 
correlation method discussed above requires simultaneous 
temperature and velocity measurements. An advantage of the 
eddy correlation method is that the kinematic regime 
associated with turbulent heat flux can be analyzed. 




( 8 ) 



Measures only the eddy diffusion coefficient for heat. 
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V. WAVE EFFECTS ON FLUX MEASUREMENT 



The wave-induced contribution to velocity signals in the 
near surface region has been shown to be much larger than 
the turbulent contribution. The effect of wave motions on 
the measurement of turbulent flux is analyzed below. 

The covariance between u and w in a region of wave 
influence, assuming turbulent and wave^induced motions are 
statistically independent, can be written: 

uw = UW + u^F. (9) 

For simplicity, u refers to only one horizontal velocity 
component. The first term is the contribution of the wave- 
induced motions and the second the contribution due to 
correlated turbulent velocity fluctuations. A similar 
expansion into wave-induced and turbulent parts can be made 
for temperature and vertical velocity covariance: 

0w = 0 W + 9 ' w' , (10) 

Since it has been shown that linear wave theory closely 
approximates observed data, U and W, as well as 0 and W 
(equation 4) are in quadrature and their contributions to 
the respective covariances should be zero. In terms of the 
cross spectrum: 

S (f) = Co (f) - iQ (f), (11) 

uw uw uw ’ 

wave motions should contribute only to the imaginary part. 



20 



Recalling the definition (_8), equation (9) can be written 



uw = Re 




o 



00 



S uw (f) df] = Re [u l w’ - i(U W + 6)] , (12) 



The contribution of the wave motions appears only in the 
(imaginary) quadrature spectrum because U and W are 90° out 
of phase. The quadrature contribution of the turbulence (6) 
will henceforth be neglected compared to UW, The simple 
overbar denotes RMS units. If, for some reason, the measured 
phase between U and W should differ from 90° a component of 
one signal will appear in phase with the other and contribute 
to the cospectrum, Such a contribution results in an error 
in the calculated flux. This error might be ignored for 
small phase differences were it not for the fact that the 
wave-induced signal is so much larger than the turbulent 
signal. 

Real phase shifts between U and W due to non-linear 
effects have been postulated. For example Yefimov and 
Khristoforov (1972) included viscosity effects in the equations 
of motion to show that a phase shift occurs at second order 
for steep waves. The phase shift attenuates rapidly with 
depth and increases with increasing frequency. Consistent 
frequency dependent phase shifts were conspicuously absent 
from even our shallowest velocity records in the entire wind 
wave range of frequencies (see for example figure 1). 

Other than unknown higher order or weak nonlinear effects, 
the only cause for the phase between horizontal and vertical 
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wave-induced velocities to differ from 90° is measurement 
error in directional response, frequency response, or improper 
instrument orientation. Non ideal directional response 
will be examined first, 

A. FLOWMETER DIRECTIONAL RESPONSE 

Expressing the true harmonic velocity component as 

V = | V | cos at (12) 

the ideal response of a vertically oriented flowmeter is 
pictured in polar and cartesian form in figure 3a. If the 
flowmeter response deviates from the cosine law, the measured 
signal will be in error by some amount e(at) which depends 
on the magnitude and phase (i,e, direction) of V. The 
measured signal becomes 

V* = | V | cosat + | V | e(at ) . (13) 

e(at) must be periodic in at and can thus be decomposed into 
Fourier components which are all harmonics of a but are not 
necessarily in phase with V. 

m 

e(at) = Yj e cos(nat + <{>). (14) 

n=l n n 

There are three forms in which the error can appear. 

'i * ^ z 

• d * ^ 

First, d> can be zero by virtue of the fact that the flow- 
meter will not respond to motions perpendicular to its axis. 
The poor directional response of the flowmeter in this case 
has the effect of broadening the bandwidth of the measured 
velocity signal by inclusion of the harmonics of e(at). 
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Figure 3a. Ideal Flowmeter Directional Response. 





Figure 3b. Non-Ideal Flowmeter Response without Phase Error. 
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Figure 3c. Non-Ideal Flowmeter Response with Phase Error, 



For flowmeters with noticeably poor directional response (as 
in figure 3b), this may create difficulties in analyzing 
velocity power spectra for slopes indicative of an inertial 
subrange above the wave frequency band. The electromagnetic 
flowmeters used here, and others tested by NOIC (1974), showed 
a typical first harmonic error response on the order of 2 to 
4 percent of the ideal response amplitude. Table I indicates 
that the magnitude of turbulent fluctuations was on the order 
of 10 percent of the wave induced amplitudes. Contamination 
of the high frequency turbulent signal by harmonics of 
e(at) might therefore be possible. 

As long as the flowmeter does not respond to perpendicular 
flow, however, the measured signal must be in phase with the 
actual motion. If, for example, the covariance is calculated 
between W* and another signal with which W should be in 
quadrature (say U), the two measured signals, including all 
harmonics of e(at), are still in quadrature. By virtue of 
their orthogonality there is no erroneous flux contribution. 
Surface waves usually occupy a frequency band of at least 
two octaves and there is a possibility that the first harmonic 
of the error signal in W* is not in quadrature with the U 
signal at that frequency. Since the phase between spectral 
wave components is random, however, the average contribution 
can be expected to be zero. 

The second case of (14) is when the <f> are zero but 

n 

e(+ ^) is not because of a D.C. bias in the system. Errors 
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of this sort are inconsequential because D,C, values are 
normally subtracted out of the fluctuating signal before 
calculations are made. 



The last schematic in figure 3 is for the case of <J> n f 0. 
This kind of response occurs when the measured signal at 



at = + ^ differs from that at at = r In other words, the 



meter response to a velocity component perpendicular to its 
axis depends on the direction of that velocity component. 

An error of this type does result in a phase shift of the 
measured signal with a consequent wave-r-induced contribution 
to the calculated covariance. Ducted type flowmeters in 
particular can respond erratically at very large angles of 
attack. Smith (1972) circumvented this difficulty by 
tilting his ducted impeller current meter array from the 
principal gravitational axis while measuring turbulent flux 
under Arctic ice. Unfortunately in the near surface region 
of the open ocean, the three dimensional wave motions have 
no preferred direction and this method does not work. 



B. FLOWMETER FREQUENCY RESPONSE 

There is an inherent phase lag in the response of any 
real sensor. Instruments for measuring turbulence are usually 
chosen to exhibit a flat amplitude response over the 
frequency band being measured. The phase response is seldom 
reported. Measurements using a variety of instruments, 
whose outputs are subsequently compared, are particularly 
susceptible to phase errors due to synchronization problems 
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and non-uniform phase response among the instruments. Later 
results will point out how extremely sensitive correlation 
measurements in the layer of wave influence are to miniscule 
errors in the measured phase. Care must be taken to ensure 
that all instruments have extremely short time constants or 
at least that the phase response of all instruments is the 
same. 

C, FLOWMETER MISALIGNMENT 

The other source of phase error in the measured velocity 
signals to be examined is that due to improper orientation 
of the velocity sensors with respect to the natural axes of 
wave motion. Since horizontal wave motions are directional, 
the orientation of the natural horizontal axes cannot be 
specified except that they are orthogonal to the local 
gravity vector. For each small band of the wave spectrum, 
however, it can be assumed that the direction is constant 
for a reasonable length of time. With this in mind, the three 
dimensional field of wave motions at a point can be represented 
by a Fourier sum of two dimensional motions, namely: 

m 

W = V |W I cos(a t + <p ) 

1 n 1 n T n 

n=l 

m 

U = y |u I sin(a n t + 4> n ) cos3 n (15) 

n=l 

m 

V = Z |U n |sin(a n t + <P n ) sing n 
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where 0 q are the angles between the horizontal wave component 
directions and some arbitrary reference direction. 

Choosing one Fourier component and setting 0 equal to 
zero allows the measurement of the motion to be simply 
analyzed by assuming that the flowmeter axes are in the plane 
of motion defined by U and W . Figure 4 shows the axes of 
the wave motion and two flowmeter component axes which are 
tilted by an arbitrary angle. From figure 4, the measured 
signals U* and W* are given by 

U* = |U I sinat cosa. . + |W j cosat cosa 01 

n 1 n 1 11 1 n 1 31 

(16) 

W* = lu I sinat cosa.„ + |W I cosat cosa 00 • 

n 1 n 1 13 1 n 1 33 

The effect of misalignment of the velocity sensors on 
turbulent flux measurements can now be calculated by taking 
the covariance between U* and W* and 0 . Recognizing that 

Fourier wave components are independent , the covariances 
U*W* and 0 W* can be summed over all wave components to give 
the total covariance contribution for the entire wave 
spectrum (still assuming all 0^ = 0). In terms of RMS 
amplitude this results in 



u*w* 


2 2 

= U cosa^ sina^g - W cosa^g sina^ 


(17) 


ew* = 


0U sinagg 


(18) 



Equations (17) and (18) are the errors due to wave motion in 
measured turbulent flux if the flowmeter axes are misaligned. 
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Figure 




4. Measurement of two dimensional wave motion by 

flowmeter in the plane of motion. If the flowmeter 



axes are orthogonal : 
a 13 90 a 33 ’ 



a 



11 



= a 



33 
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= 90 + a 
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From equations (9) and (.10), the percent error introduced 



into the flux measurement is 



u*w* - u f w* U*W* U^cos .^sing^-W^cos ggSin 
u' w ' u T w ' u * w v 

0w* - 0'w’ 0U sina^g 

0 'w' 9 ' w* 



(19) 



( 20 ) 



The true value of the turbulent covariance is not known, but 
the turbulence correlation coefficient (r*) can be defined: 



r^ w = u'w7(u' w l ) < 1 (21) 

r 0 W = 0 vce -1 w l ) < 1 (22) 

Dividing numerator and denominator of (19) and (20) by u' w' 
or 0 ' w' as appropriate gives the percent error in terms of 
the tilt angles and the turbulence correlation coefficient. 
Equation (20), for example becomes 

0 U sina™ 

Percent error = — — — . (23) 

0’w' r’ 

0w 

Equation (23) states that the percent error is inversely 
proportional to the "turbulent signal" to "wave-induced noise" 
ratio and the correlation between turbulent components. It 
is directly proportional to the tilt angle. Figure 5 shows 
the percent error (multiplied by the signal to noise ratio) 
for various values of the tilt angle and turbulence 
correlation coefficient. 
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Figure 5. Turbulent Flux Error Due to Misalignment of Flowmeters, 



Equation (.19) can he simplified for illustrative purposes 
by assuming that the flowmeter axes are orthogonal so that 
ot 11 = a 33 - Further assuming that the tilt angles are small, 
such that cosa^ 1, equation (19) becomes 



U*W* CU 2 - W 2 ) sina 33 
u‘w' uV 

and the expression for percent error can be written 



(24) 



CU 2 r K 2 ) sina 33 
Percent error = 



(25) 



u'w' r' 

uw 

Figure 5 is again applicable if the turbulent signal to 
wave-induced noise ratio is interpreted as 



u 'w* 




Equation (25) shows that, for the restricted case described 
above, the error in deep water (where W = U) should be zero. 

An example of the momentum flux error can be calculated 
using the measured intensities of turbulent and wave-induced 
motions from Table 1. Assume a tilt angle of 5° and a 
turbulence correlation coefficient, r' , of 0.2, which has 
been borrowed from Bowden and Howe's (1963) results in an 
estuary. The turbulent signal to wave-induced noise ratio 
for run 22 is obtained from Table 1 and is equal to 0.064, 
Entering figure 5, the product of the percent error and 
signal to noise ratio is 44. Solving for the percent error 
gives a value of 688%. 
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Since there are no representative values available for 
the correlation coefficient between turbulent temperature and 
vertical velocity fluctuations, estimates of the heat flux 
error rest on the assumption that r^ w is approximately equal 
to Using the previous values of 0.2 for r* w> 5° for the 

tilt angle, and the signal to noise ratio from Tables 1 and 
2 CO, 29) results in a percent error for the heat flux measure- 
ment of 155%, The heat flux measurement is less sensitive 
to tilt errors because of the larger signal to noise ratio. 

The accuracy with which flowmeters must be aligned to 
achieve an acceptable percent error can also be calculated. 

For a 100% error, a signal to noise ratio of 0,064 and an 
assumed turbulence correlation coefficient of 0,2 the 
required alignment accuracy for momentum flux measurements 
is + 0,73°, For the heat flux measurement, using the same 
values as the previous example, the required alignment 
accuracy is 3.3°, 

For the general case where the horizontal axis of the 
flowmeter does not lie in the plane of all orbital motions 
(due to the directionality of the waves), the above develop- 
ment becomes more complicated. In this case, let the 
reference direction be the azimuth of the horizontal flowmeter 
axis so that 3 n are the angles the horizontal wave motions 
make with the flowmeter axis, as in figure 6. With the 
coordinates defined as above all wave motions are in the 
1-3 plane and V is zero for each individual spectral component. 
The response of the flowmeter to a typical wave field of 
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Figure 6. Measurement of Directional Two Dimensional 
Wave Motion by a Flowmeter not in the Plane 
of motion. 
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random direction is 



(U* V* W*) = (U W ) 
n n n' v n n y 



c° s COii) n cos(a 12 ) n cos(a 13 ) n 



CosC °31 ) n COS(a 32 ) n COs(a 33 ) n / 



(28) 



where the last term on the right consists of the directional 

cosines between the unprimed and primed axes. The angles 

(a^ ) n now include the effect of wave direction and may be 

different for different Fourier components. Considering a 

single wave component making an angle 8 Q .with the azimuth of 

the horizontal flowmeter axis, the covariances between U* , V* 

’ n n 

and W* and between 0 and W* are easily calculated 
n n n 



U*W* = U cosa^cosa., W cosa 01 cosa 0 „ 
n n n 11 13 n 31 33 



V*W* = -U cosa,„cosa 10 + \V cosa„ 0 cosa„„ 
n n n 12 lo n 66 66 



(29) 

(30) 



0 W* = 0 U cosa 00 
n n n n 33 



The summation of wave components can not be accomplished, 
however, without knowledge of the direction spectrum of the 
waves. Even then it is difficult. 

Again assuming that the flowmeter axes are orthogonal, 
some statements can be made about the wave-induced error. 

It is clear that the error contribution to the U* signal is 



w n sln(a 33>n 

regardless of the direction of the waves. On the other hand, 

the error contribution to W* is a maximum at 8 = 0 (i.e. 

n n 
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when U q is largest) and decreases as the quantity cos$ n , 

The previous result for deep water thus does not hold and 
some error is to be expected in deep water flux measurements. 
The error for B n = 0 in deep water is zero, as in equation 
(24), But the error for (3^ = 90° is a maximum given by 

- c 0 sCc. 33 ) n sin( 033 ) n , (31) 

since the U* signal consists only of the projection onto 
the 1* axis. 



VI, DESCRIPTION OF THE EXPERIMENT 

The Naval Postgraduate School CNPS) has conducted a 
series of turbulence measurements using a fixed array of 
sensors mounted on the NUC oceanographic research tower off 
San Diego. The array was designed to simultaneously measure 
surface elevation temperature, water particle velocity, 
conductivity and sound transmission parameters in a two cubic 
meter volume of ocean. A detailed description of the sensors 
and ancillary electronics was given by Hagen (1974), The 
data reported here consists of the output of a resistance 
wire wave gauge, a two-component electromagnetic flowmeter, 
and a thermistor adjacent to the flowmeter. The instruments 
were all mounted in a vertical line. The vertical separation 
between thermistor and flowmeter was six centimeters. Analog 
sensor outputs were recorded on FM multichannel tape 
recorders and subsequently passed through matched analog 
lowpass filters before being digitized. Prior to further 
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analysis, all records were highpass filtered using a running 
average filter with a 0.00675 Hz cutoff frequency to subtract 
out the low frequency fluctuations due to internal waves. 

Data were obtained in water approximately twenty meters 
deep. A persistent thermocline was observed at mid-depth, 
the mean temperature profile being characterized by two 
nearly isothermal layers separated by a relatively sharp 
gradient. Figure 7 shows the position of the sensors in 
the water column and relative to the thermal structure as 
sampled by bathythermograph (Hagen, 1974). The thermocline 
was found to correspond closely with the 12° isotherm in 
most cases. Internal wave activity was observed during each 
run . 

The thermocline, when internal waves are present, is a 
region of relatively high velocity shear combined with the 
large temperature gradient. Velocity shear due to internal 
waves could not be determined with any certainty because the 
internal wave velocities were only slightly higher than the 
noise level of the flowmeters. The maximum apparent shear 
was on the order of 1 cm sec~^m~^. The uncertainties in 
low velocity water motion measurements precluded estimation 
of the Richardson number. 

Turbulent heat fluxes were calculated for six runs made 
at various depths. The array was positioned in each case 
such that the sensors were near and, in some cases, in the 
thermocline. As internal waves passed the array, the 
Lagrangian flow field was sampled to a vertical extent of 
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LOCAL TIME 



about 2^3 meters. Thus, the averages calculated from the 
experimental records were spatial as well as time averages. 

The device of filtering out low frequency temperature 
and velocity data merely prevented internal wave components 
from biasing the measurements. Since available averaging 
times were insufficient to provide proper estimates at 
internal wave frequencies, it was decided to ignore them. 
Internal wave effects, however, were still observable in the 
filtered records because of the above mentioned sampling 
of the Lagrangian temperature and velocity fields. 

The calculated turbulent heat fluxes were somewhat 
higher than expected. Advection of heat away from the 
experiment site was discounted as a reason for this result 
due to the low current velocities in the area. A counter- 
balancing flux contribution at the internal wave frequencies 
which had been filtered out was considered. The flux was 
calculated from unfiltered records but nothing conclusive 
could be observed due to lack of sufficient record length. 
The imposed vertical scale of twenty meters and the stable 
density gradient were presumed to limit turbulence to 
frequencies higher than that of the internal waves 
(Kitaigorodskii , 1973). Attention was subsequently focused 
on the extremely large effect of the surface waves on the 
temperature and velocity records. 
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VII, MEASURED HEAT FLUX 



The heat flux measurements made from the NPS array were 
suspected of containing significant errors due to flowmeter 
misalignment. In light of the preceding error analysis it 
was considered possible to make certain corrections to the 
flux measurements if the tilt angle of the vertical velocity 
sensor were known. The cross spectrum between surface dis- 
placement and vertical velocity was available, and with the 
appropriate approximations, a relationship between the 
measured phase angle and the actual tilt angle was derived. 

Several fortunate circumstances combined to allow the 
tilt angle to be calculated. The physical mounting of the 
flowmeters was such that rotation was possible only in the 
l'-3' plane (.ref erring to figure 6). In addition, coastal 
environments are usually characterized by relatively narrow 
banded wave direction spectra so the assumption could be made 
that the wave motions were essentially in one direction. 

Since the predominant wave direction was consistently within 
30° of the azimuth of the plane of the l'-3 ? flowmeter 
components, small angle approximations could be used. Refer- 
ring again to figure 6, for small can be considered 

equal to its projection on the 1-3 plane and the problem is 
simplified to that illustrated in figure 4, where 

a 33 90 " a 13‘ 
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The cross spectrum between wave displacement and vertical 
velocity for this situation is derived in Appendix B and 
can be written 



2 




S nw Cf) df = 2 



CCo(f) 



i Q( f ) ) df = 



= Ch k w' + n U sina 33 ) ^ i(jl W cosa 33 ) 



(32) 



Making the assumption that random surface displacements 
(e,g, breaking wave crests) are essentially uncorrelated 
with turbulent velocities at a depth of 5 meters (the minimum 
measurement depth), or greater, allows the h’w' term to be 
neglected. The vertical integral scale, measured from 
temperature data to test this assumption, was 1,5 meters. 

The cross spectrum is now defined totally by wave-induced 
motions and the phase angle in the wave band is 



<j> = tan 



-1 




Q(f ) df 



Co( f ) df 



T1 W cosa 



n U sina 



33 

33 



% 



Solving for a 



33’ 



a 33 = tan 



-1 



W 

U tan <f> 



(33) 



(34) 
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The cross spectra between waves and vertical velocity were 
obtained for the two shallowest runs and was calculated. 
The results are given in Table III, 

Run 6 : seas from 310Tj swell from 290T 
Run 10: seas from 30QT; swell from 270T 
vertical flowmeter Run 6 Run 10 

tilt angle (.<* 33 ) 3,54° 3,62° 

TABLE III, Tilt angle of vertical flowmeter axis. The 

azimuth of the plane of vertical and horizontal 
flowmeter components was 270T, 



Using a value of equal to 3.6 for all runs, equation 

(18) was used to calculate the wave-induced error and the 
correction was applied to the measured heat flux. The 
average temperature gradients during each run were also 
measured from thermistors approximately one half meter above 
and below the point of flux measurement. The eddy diffusion 
coefficient was then calculated from 



k h - 



0 ’w' 



39/3, 



(35) 



The turbulence correlation coefficient (r’) was calculated 

v ew 



from corrected values of 0 f w’ and the turbulent intensities 
from Tables I and II. The results are given in Table IV. 

The turbulence correlation coefficients were all less than 
one except for run 26. A turbulence correlation coefficient 
greater than one indicates that there were errors in the 
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measurement which were not corrected for and, hence, run 26 
is disregarded in the discussion below. 

The variations of heat flux as internal waves pass a 
fixed sensor should reflect the variations in velocity and 
temperature gradients encountered. In the absence of large 
scale mixing events, such as the breaking of internal wave 
crests, the heat flux should be proportional to the temper- 
ature and velocity gradients. It was therefore expected 
that the measured heat flux would increase as the thermocline 
region was more extensively sampled. Figure 7 and Table IV 
show that this expectation was not fully realized, since a 
simple relationship between temperature gradient and heat 
flux was not observed. 

Of the runs conducted at night (17, 18, 19), the first 
two indicated an upward transport of heat presumably due to 
convective overturn. Both these runs, although close to the 
thermocline, mainly sampled the isothermal regions above and 
below it. It was not immediately apparent why run 18, below 
the thermocline, would show an upward heat flux, because 
surface water would have had to be cooled quite radically to 
sink through the thermocline. The change in shape of the 
temperature profile before and after the run indicated that 
a large scale mixing event might have occured during the run, 
which would account for the upward heat flux. Run 19 sampled 
a region of relatively high temperature gradient which 
resulted in a net downward flux, as expected. The heat flux 
was less, however, than for a daytime run in a region of 
similar temperature gradient (run 28a). 
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Runs 22 and 28 were conducted during the late morning 
and afternoon, respectively. Run 22 sampled the thermocline 
extensively and there was evidence of a large scale mixing 
event, both in the temperature time series and in the BT 
samples. The heat flux was consequently large. Run 28 
was divided into two runs because the first part of the run 
sampled an isothermal region while two very distinct internal 
wave troughs passed by the sensors during the second part of 
the run. There was no evidence of a large scale mixing 
event occuring at the level of the flux measuring sensors 
in run 28, although the thermistor placed a meter below did 
record a possible mixing event. This run was considered to 
be measuring only small scale turbulent flux at various 
distances from the thermocline. 

The high turbulence correlation coefficients observed for 
runs 19, 22, and 28a correspond to relatively long periods 
of time during which the sensors were in the thermocline 
region. The inference can be made that gradient type eddy 
diffusion in this region is an efficient heat transport 
mechanism which can probably be closely approximated by a 
gradient diffusion model of the form of equation (35). The 
relationship between K„ and other thermocline parameters 
(30/3z, 3U/3z) could not be determined because only two of 
the six runs (19, 28a) did not show evidence of larger scale 
mixing events. 

On a qualitative basis, the results showed that small 
scale turbulent heat transport in the upper ocean is not 
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easily modeled. The mechanism for transfer seems to be a 
combination of gradient type diffusion and convection. Mixing 
at large scales is a factor and may be due to breaking of 
internal wave crests. It appears that large scale processes 
reestablish sharp temperature gradients which are again 
eroded by eddy diffusion and convection. Large increases 
in turbulent kinetic energy were not observed in the thermo- 
Cline region, indicating that production of turbulence by 
internal wave-induced shear was not as important as produc- 
tion at the surface. The local Richardson number must there- 
fore have been greater than the critical value of 0.25 for 
the maintenance of shear-induced turbulence, 

VIII, CONCLUSIONS 

Analysis of near surface velocity records in shallow 
water indicated that the ratio of turbulent to wave-induced 
kinetic energy was of order 10^^. Because of the relatively 
large wave-induced fluctuations, large turbulent flux 
measurement errors can occur for small errors in directional 
response, phase response, or instrument alignment. Alignment 
errors are by far the most troublesome because they cannot 
be determined in the laboratory before measurements are made. 
Flowmeter misalignment of 3,6° caused errors in the calculated 
heat flux of about 25% to 100%. Momentum flux measurements 
would have been affected to a greater degree. It was shown 
how momentum flux measurements in the open ocean would be 
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less sensitive to misalignment errors than shallow water 
measurements. Deep water misalignment errors would be 
minimized by aligning the flowmeter axes as near as possible 
to the predominant direction of wave motion. 

The heat flux measurements reported here were only an 
initial effort and the results could certainly be improved 
upon. Verification of the applicability of the gradient eddy 
diffusion model in the thermocline was not possible from the 
sparse data. It is felt that analysis of the remaining NUC 
tower runs will provide sufficient data to establish a 
relationship between the temperature and velocity structure 
in the thermocline and the eddy diffusivity. 
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APPENDIX A 



The relative magnitudes of the wave-induced and turbulent 
fluctuations can be approximated by assuming that: 

1. ) the wave-induced and turbulent fluctuations are 

independent , and 

2, ) the wave-induced fluctuations can be described by 

linear wave theory. 

Using the first assumption, the spectra of (1) 3-re related 
according to: 

S u Cf) = S 0 (f) + S u ,Cf) (a) 



The second assumption allows S^Cf) to be written in terms of 
the spectrum of the surface displacement (S^(f)) as: 

SyCf) = | H 2 ( f ) | S n (f) (b) 



where H(f) is a transfer function derived from linear theory. 
The coherence between U and n is given by 



r Un (f } 



s un (f) 



SuW 



s n (f) 



= l 



(o) 



where the identity holds for a linear constant parameter 
system. Now 

|s Un ( f )| 2 = SjjCf) S n (f) (d) 

and by assumption (1.) 

I s un <1)|2 " |s un (f>|2 ” s v (1) s n (f) <e) 
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The measured coherence between u and n is 

2 I s un ( f ) I 2 

(f) 1 

1 uri 



s u (f) s n (f) 



substituting (e) and (a), 



2 S_Cf) S (f) - S ,(f) 

x (f) = -a = -a a 

ur) ' 



S Cf) 

u v * 



S (f) 
u 



which results in 



S ,(f) = S (f) (1 - y ) , 

u » v / u v y v 1 UTl 



(f) 



(g) 
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APPENDIX B 



Cross Spectrum Between Fourier Components of n and W* 



f - 7 



*u,> i 



W*= LT^Un o-± cot o< + Wcoto-i ccxjcv 

= i U A+n (cr-t ) + £ U ^in (<rk -o < /3 ) 

~ ^ CCK3 (vk - c<j 3 ) + "W cos (cri. t ) 

^■}JW+ ('?) y. ^ £ -k y IT Ain (o-t +crv 

* i 7 W^incri ^Un (cr-t - ck / s> ) 

-c TV Ac/: cri cos (cri -r <jZ - ^ 33 ) 

+ & cri Cot (cri: + cT 

• diz 

ftyv* (?) ~ A J TT cosck J3 axs (~az) 

** £ ^ TV cot c<jj C-<tt) 

-S' -ijw* T J ~ 

\ 

I / ^ ? x 

^ J {iyUco^ oc ;J [«* err ^ <3ttJt - l cos ctTaU cSrrjr\ 

' C ; 

- it y Tv COI [ sU*t(j'C CCS Att ^TT - L A^urf ALonATTfr]\ 



'olT 






S: 



)V* <f ~ 0 7> / ^ ^^CT 



^ v- *• ~j/: "J TaJ COS <\* j j 

othca cut sc. . 
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Cross Spectrum Between Fourie 



r Components of U* and W* . 



U — LT -nTt crt Conor,, *■ COO o’-i COn <x 3l 

W* - w < j-£ Codoc, 3 ** W CGQG-L CtX ’ 0^33 

U* = £ U xTifat +«,,) + ±17 ±in(o~t -«,. ) 

* * Wc ^(vt r i-Wcon(crt 

W* - ~j[ U rdnCo'i + o(, 3 ) + a If Ccr-t -«,j ) 

** •? W coo (cri -©<33 ) ^ ^ ZV con (&~t /or J3 J 

^Q'U* ^ ~ 7 " fl ('crt r-cn,,) nfn Co-t-'o-T+^jJ 

/1 

/ ^ U ( rt +«„) ^ (at -+o-T ~<K 3 ) 
u CTl V. Tn (crt -nx,, ) con (crt * cyt - °< iS ) 

~f Till 4on (<Jt + x„) con (crt ^OZ-t-Xa) 

+ -£ U ('erf -<*„) ^ on (°^ r crT ^ *</ j J 

+ ~4 1?^ (crt -o<„ ) Mn (crt •*■ (Ttr -oc ti ) 

* ~J? IW xUn (crt -<*„ ^ cr« / err -** 33 ^ 

•/ TJtJ ^fn (o"i-x „ ) cor) (erb + ot ~rCX i3 ) 

/ -y - R/77 coo (cri -(><£ ) (cr-£ + crt + 

v tiru coo (<yt ~oc 3 ,) (crt /err -ck /3 ) 

•+ ~y W x con (crt ~o< 3l J con (crt +<rX -<x j3 ) 

-*■ "y W 2 coo (crt -Xu) con ('crt +C-T ■* * 33 ) 

-* -y- E7C7 coo (o-f ^ (erf /err 

~t TV T con (cri oo(jf ) (art rerv - o< /3 ) 

■*■ y- coo (ferf /<Vj, J ^ (crt + cr r -Kjj J 
/ y' £ 1/ 1 Cen (cr£ • ofj, ^ coo (<r t v-cr r / o^j ) J Tf 
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(r) - 



U ^ <rV «. -C U*u& («,, -<*, 3 ) C 04 err 

U'AJ cck (<x u +c < 33 ) Aajt\ ert - UVJ ccx ((x lt - o( 3t ) ajsa o-r 
UW cxn +« l3 ) M*\&r +4 ZTW co-i (o<ji -c*, 3 ) MaxctX 

S U ' V »(p — 7 U W C®<„ +«^, 3 ) ■*■ -f U*coi (c* H ~cx^) 

■+ T V 1 G&l ( CKji + c^ 33 )+ ~t TV *003 ((Xj- 0< 33 ) 

•f i £ gT t/ W CAJ (V/, + C^ 33 ) + -f-UW CCnj ( <V„ -0^3 J j 

# j *v" o« {exj, + <^ l3 ) - y‘U/v Cxx /o^j, - c^ l3 ) J 
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